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Abstract 

Inspired by the decaying dark matter (DM) which can explain cosmic ray anomalies naturally, 
we consider the supersymmetric Standard Model with three right-handed neutrinos (RHNs) and 
/2-parity, and introduce a TeV-scale DM sector with two fields 01,2 and a Z% discrete symmetry. 
The DM sector only interacts with the RHNs via a very heavy field exchange and then we can 
explain the cosmic ray anomalies. With the second right-handed neutrino N2 dominant seesaw 
mechanism at the low scale around 10 4 GeV, we show that (f>i^ can obtain the vacuum expectation 
values around the TeV scale, and then the lightest state from <f>\ 2 is the decay DM with lifetime 
around ~ 10 26 s. In particular, the DM very long lifetime is related to the tiny neutrino masses, 
and the dominant DM decay channels to u and r are related to the approximate u — r symmetry. 
Furthermore, the correct DM relic density can be obtained via the freeze-in mechanism, the small- 
scale problem for power spectrum can be solved due to the decays of the -R-parity odd meta-stable 
states in the DM sector, and the baryon asymmetry can be generated via the soft leptogensis. 

PACS numbers: 12.60.Jv, 14.70.Pw, 95.35.+d 
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I. INTRODUCTION AND MOTIVATION 

The cosmic ray anomalies observed by PAMELA and Fermi-LAT [l, 2} strongly indicated 
that the dark matter (DM) particles annihilate or decay dominantly into the leptons. Al- 
though the large annihilation cross sections can be realized via the Sommerfield enhancement 
or Breit-Wigner mechanism jsj, the HESS obervation 
rays gives strong constraints on the annihilation DM scenario 
with a lifetime at the order (9(10 26 )s is another elegant way to explain the cosmic ray anoma- 
lies. In particular, the contraints from the Galactic center gamma rays are much weaker jf| . 
However, the ultimate long lifetime of decaying DM becomes a non-trivial problem since the 
symmetry, which makes the DM stable, must be broken tinily. 

Supersymmetry naturally solve the gauge hiearchy problem in the Standard Model (SM). 
Gauge coupling unification in the Minimal Supersymmetric Standard Model (MSSM) implies 
the Grand Unified Theories (GUTs) at the GUT scale M GUT around 2 x 10 16 GeV. Thus, 
the DM decays via the dimension-six operators suppressed by the GUT scale is a rather 
simple solution to the long lifetime of decaying DM, and it may provide a way to probe 
the GUT scale physics sjj. Another problem in the decaying DM is how to understand the 
DM dominant leptonic decays, especially to the /1 and r final states. Without automatically 
kinematical suppressions like the annihilation models [sj], one has to employ some special 
symmetry so that the DM interacts strongly with the second or third families of the charged 
leptons, for instance, flavor Froggat-Nielson symmetry [lo| . 

Because of the DM leptonic decays, one may conjecture that the DM sector only interacts 
with the lepton sector 11] . Note that the neutrino masses are very tiny, we can parametrize 
the small couplings for the DM decay as the ratio between the light neutrino mass m v ~ 10~ n 
GeV and the GUT scale 

A ~ ~ ID" 27 . (1) 

Interestingly, this is the typical order of tiny coupling parameter rendering the lifetime 
around 10 26 s for a TeV-scale DM. Thus, it also implies the deep connection between the 
decaying DM with long lifetime and the active neutrinos with tiny masses. 

In this paper, we consider the supersymmetric Standard Models with i?-parity (R p ) and 
three right-handed neutrinos (RHNs) N{ where the neutrino masses are generated via the 
low-scale seesaw mechanism. In the DM sector, we introduce the SM singlet DM fields <pi^ 
and a discrete Z3 symmetry under which the term 0102 is invariant. At the leading order 
they couple to leptonic sector through a dimension-seven operator 

x x 4 x * i ^te (L ' A,(L ^" ) )' (2) 



2 



where coefficient C-j come from neutrino Dirac Yukawa couplings, constrained by light neu- 
trino masses in seesaw mechanism. It can be obtained after integrating out the heavy 
right-handed neutrinos (RHNs) and a superheavy field X with mass A, provided that fa 
only interacts with the RHNs mediated by X. Note that there is a GUT scale in GUTs, we 
shall assume A ~ Mqvt- Interestingly, this is the exact scale that we need to produce the 
correct DM density via the freeze-in mechanism. In our models, the vacuum expectation 
values (VEVs) for </> 1>2 are close to the RHN masses. Thus, Eq. ([2]) gives the small parameter 
A approximately if we identify the terms in the bracket as m u . Furthermore, the preferred 
fi and r decay channels can be related to the neutrino tri-bimaximal mixin g (T BM) sce- 
nario with the second right-handed neutrino N2 dominant seesaw mechanism 12j|. In short, 
cosmic ray anomalies, if confirmed further, potentially have deep correlation with neutrino 
physics, especially such seesaw mechanism. 

The DM relic density generally conflicts with standard thermal freeze-out scenario in the 
decaying DM models if the scalar components of the DM fields acquire VEVs. It is not diffi- 
cult to understand it from the effective operators 0i</>2-£ 2 (£ are the operators for particles in 
the MSSM) which can generate large annihilating rate. However, they also catastrophically 
make fa decay very fast. So one usually considered the non-thermal DM production, for 
example, a detail study given in Ref. 13j. Unlike the weakly interacting massive particle 
(WIMP) scenario, the non-thermal production usually loses the natural predication on DM 
abundance. Recently, it was proposed that the feebly interacting massive particle (FIMP) 
may be an alternative to WIMP [l^], shedding light on decaying DM. Typically, the FIMP 
involving a coupling at the strength (9(10~ 13 ) ~ TeV / 'Mqxjt for the decay dominated freeze- 
in mechanism, or a larger one ~ O(10~ 11 ) for scattering dominated freeze-in mechanism. 
Amazingly, in our model, fap. must couple to the RHNs by the dimension-five operators 
suppressed by Mx somewhat smaller than Mqut and by ~ 10 4 GeV, based on proper 
decaying lifetime of DM. Similar results hold for the scattering process dominated freeze-in 
mechanism. Therefore, in our decaying DM model, its relic density again is a "miracle" via 
the freeze-in mechanism. 

As a by product in supersymmetric SMs with freeze-in mechanism, we are able to solve 
the small scale problem for power spectrum, in the presence of a metastable R p — odd state 4> 
in the DM sector. The point is the following: the whole supermultiplets fa^ are freezed into 
the thermal bath. We assume that vn^ R + m-r, > m^, where <pn is the lightest state and the 
DM particle while 0' is the lighter R p — odd state and has a mass close to the DM particle 
<Pr. By virtue of x R p , the leading decay mode of (ft is to the lightest supersymmetric 
particle (LSP) plus <pR. Provided that and the LSP are respectively sufficient heavy 
and light, the relativistic LSP is produced from late decay at some sufficient late time 
17 ~ 10s — 1000s. This warm DM component is just the key to reduce power spectrum on 
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small scale [15]. 



In addition, we can still explain the baryon asymmetry via soft leptogenesis 16| since 
the seesaw scale is low around 10 4 GeV. In supersymmetric seesaw framework, when the 
new CP- violating sources in the soft terms dominate the sneutrino(s) N CP- violating decay, 
the so-called soft leptogensis 0, 0| is indeed able Jo produce enough lepton numbers in 



our model. Before the discovery of gaugino effect [32|, soft leptogensis suffers the highly 
suppressed bilinear soft mass term B N .M Ni NiNi with B N . < 10~ 3 M S usy where M S usy is 
the universal supersymmetry breaking scale. In this paper, we assume that the trilinear soft 
terms AY Nl3 ' NiLjH u are the only sources of CP-violation in the supersymmetry breaking 
soft terms. Interestingly, enough baryon number density can be generated naturally. Even 
the baryon number density tends to be overproduced if M^ 2 is too light ~ O (TeV), we can 
choose a relatively smaller universal A term or tune its CP violation phase to obtain the 
observed baryon asymmetry. 

This paper is organized as follows. In Section II, we present the model, and discuss the 
decay of DM, its relic density, as well as the relation between neutrino physics and cosmic 
ray anomalies. In Section III, we study the phenomenological consequences of our model, 
such as the solution to the small scale structure problem, and the low-scale soft leptogensis. 
In the Appendix A, we briefly review the freeze-in mechanism. 

II. THE SUPERSYMMETRIC DECAYING DARK MATTER MODEL WITH N 2 
DOMINANT SEESAW MECHANISM 

A. The Decaying Dark Matter Model 

We consider the supersymmetric SM with three RHNs and -R-parity, and introduce a DM 
sector. As we know, a well defined DM sector should not only have a DM particle at the TeV 
scale, but also spontaneously breaks the discrete symmetry that stabilize the DM particle. 
The simplest dark sector contains a SM singlet and a Z 2 discrete symmetry under which 
only is odd. To have a decaying DM, we break the Z 2 symmetry by giving a VEV to 
0, i.e., ((f)) ^ 0. Thus, at the leading order, couples to the observable sector through 
dimension-five operators 2 A r J 2 /A, which can be obtained from a renormalizable theory after 
integrating out a heavy field with mass A. This decaying DM can explain the cosmic ray 
anomalies and satisfy the other phenomenological requirements that we shall consider. For 



an example, see Ref. 19j. However, in such an simple model it is very difficult to break 
Z 2 symmetry naturally. Concretely speaking, we may have to introduce another SM singlet 
field S' that couples to via a superpotential term S'(f> 2 . This superpotential term provides 
the quartic term to the scalar potential of 0, and then we can realize the Z 2 symmetry 
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breaking. After Z 2 symmetry breaking, and 5" will mix with each other. Because we 
cannot forbid the direct couplings between the SM singlet S' and the observable sector, this 
simplest model is excluded unless we have huge fine-tuning. 

Therefore, we consider a DM sector with two SM singlet fields 0^2 and a discrete Z 3 
symmetry. Under Z3 symmetry, 0^2 transform as follows 

01 -> W0i , 02 -> W 2 (p2 , (3) 

where w = e l2n ^ 3 . All the other fields in our model are neutral under the Z3 symmetry, so, 
any renormalizable coupling term between the DM sector and observable sector is forbidden 
due to the Z3 symmetry and .R-parity. In particular, the traditional particle physics in the 
observable sector will not be affected. The most general Z 3 — invariant and renormalizable 
superpotential in the DM sector, as well as the corresponding supersymmetry breaking soft 
terms are 

W DM =y0? + y0^ + M 0!02 , 

^JJ^I 2 + mj 2 |0 2 | 2 + (%A X ^ + ^A 2 0| + ^ikr 0!0 2 + /i.c.^) . (4) 

In fact, this model not only preserves Z%— symmetry, but also a trivial -R-parity. In the 
following, we shall prove that this simple DM sector can break the Z3— symmetry sponta- 
neously, and has a proper spectrum with a TeV-scale decaying DM coupling to the observable 
sector. 

As poingted out in the Introduction, to have the desirable DM lifetime, abundance and 
decay products, the DM should couple to the RHNs via the dimension-five operators sup- 
pressed by the GUT scale ~ Mqut- This can be achieved by integrating out a heavy SM 
singlet field X, which mediates the interactions between the DM sector and observable 
sector. So we consider the following superpotential 

W D^N? + YgL&Nj + ^fXN.N, 

( M x \ 
+ Ax</>0i02^ + ( — ^— A + irrevelant terms 1 , (5) 

with Mx ~ Mqut- For simplicity, we have assumed that the RHNs are in the mass basis. 
We can explain the neutrino masses and mixings by employing some non-Abelian flavor 



symmetry such as A4 [20J, although we do not consider it here. In addition, we do not 
consider the superpotential XH u Hd so that we can explain the PAMELA experiment. This 
can be realized in the five- dimensional scenario compactified on S l jZi (or in the M-theory on 
S l /Z 2 ) where X and Hu/H^ are localized on the different D3-branes on the two boundaries 
of S 1 /Z2 while the right-handed neutrinos are in the bulk. 
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To obtain the effective action below the scale M x , we integrate out the heavy field X 
through its equation of motion 

M X X + \x^i<h + *XiNf = . (6) 

So we obtain the desirable dimension-five operators, which describe the interactions between 
the DM sector and RHNs. The effective superpotential are 

W N , eff =W hld - 0102^ 

+ (M^ N 2 + Y » Li H u N 3 >J + (...), (7) 

where dots denote the irrelevant corrections after integrating out X. Also, the corresponding 
supersymmetry breaking soft terms are given by 

'hid 



-C soft DC^ + mjfjNi 



P'P ^ ~ R ^ 

+ ( -z-AvrtthfaNi + AjYifLiHuNj + -£-M Ni N? + h.c. ) , (8) 



where 



r <t> _ ^x^xi f q\ 

Especially, the effective action described by Eqs. ([7]) and (jSJ) will provide the dynamics to 
freeze-in the DM particles, and generate baryon asymmetry at the scale around M^.. 

To obtain the effective action below the right-handed neutrino mass scale, we further 
integrate out the RHNs through their equations of motion 

M Nj N, + YgLiH* - ^xifafaNi = . (10) 

Then, the leading order operators coupling 0j to the SM particles are dimension- seven 
operators presented in the Introduction, Cij(f>i(f>2(LiH u )(LjH u ) (also see Fig. [Q, which can 
account for cosmic ray anomalies. The operator coefficients are 

\ \ Y N Y N 
ij ~ M x M\ ' 

where I is summed over three RHNs. Those coefficients have a clear correlation with neutrino 
Dirac Yukawa couplings as well as the light neutrino Majorama mass matrix. After Z3 
symmetry breaking by the VEVs of 0j, we obtain the renormalizable interactions between the 
DM and (s)leptons from superpotential terms <fiviVj, as well as the dimension- five interactions 
4>Vi(LjH u ). The dimension- five operators are interesting since the DM particles can decay 
dominantly to the /i and r leptons due to the neutrino TBM. 
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FIG. 1: Feymann diagram for the dimension-7 operators Cij(j}i<p2(LiH u )(LjH u ) /Mx generated 
by integrating out X and iVj at tree level. 

B. Spontaneously Z3 Symmetry Breaking and Decaying DM 

Cosmic ray anomalies can be explained elegantly by the long-lived decaying DM with 
lifetime r ~ 10 26 s that decay dominantly to the charged leptons. Because the DM lifetime 
is so long, it is natural to have a symmetry if the DM is stable. In our model, this symmetry 
is the discrete Z3 symmetry. To break the Z3 symmetry spontaneously, we consider the 
relevant scalar potnetial V(<f>i) from Eq. (j3J) 

V{<t>i) =|Ai0? + M 2 | 2 + |A 2 2 + M^l 2 

+ mJJ&l 3 + mjj0 2 | 2 + + ^A a $ + B+M+<hh + h.c^ . (12) 

Note that Mx 3> M^, the contributions to the low energy effective scalar potential from 
the superpotential X0!0 2 are very small, and then we do not consider them. Because 
Eq. (j!2p contains quite a few parameters, analytical study is pretty difficult. To reduce the 
parameters in the DM sector, we assume that the squared soft masses m? are universal, 
and the trilinear soft terms A\ i are universal. Moreover, to avoid the Landau pole problem 
for Yukawa couplings below the GUT scale, we choose Ai = A 2 = 0.3. 
First, we parametrize the fields 0; as follows 

— 7f — ' 

where "0" denotes the interaction eigenstates. We require that the spectrum have the 
following properties: (i) The lightest scalar as the DM particle should be about 2 TeV from 
the Fermi-LAT electron excess at high energy region, (ii) There should be a heavy and 
sufficient long lived R p — odd fermion with mass about 5 TeV so that we can solve the small 
scale problem. Although these requirements impose some constraints on parameter space, 
they can still be satisfied easily. In the following, we present an explicit example whose 
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input parameters are 



A 1 = A 2 = 0.3, 1.0 TeV, 

mj. = 200 GeV 2 , = -600 GeV, A x . = 600 GeV. (14) 

Let us comment on the above choice of parameters. With the fixed Aj, without tuning 
on supersymmetry breaking soft terms, larger will generate larger VEVs for fa as well 
as heavier spectrum, which is disfavored by the DM mass requirement. Because A\ { have 
the same sign, we choose negative so that the VEVs for fa have the same sign as well 
and we can have an absolute stable vacuum. The larger A— terms help us to have a more 
phenomenologically interesting spectrum, i.e., to increase the mass splitting between and 
4>r and meanwhile to keep < m ( j >R + mr,. 
Numerically, one global minimum is located at 

(fa) = vi « -6.06 TeV, (fa) = v 2 » -6.57 TeV. (15) 

At this vacuum, the various mass eigenvalues and corresponding eigenstates are respectively 
given by 

m^ a « 2.60 TeV, « 4.81 TeV, 

rrifr » 2.91 TeV, « 4.83 TeV, 

« -4.80 TeV, m^, « -2.78 TeV. (16) 
<p R = 0.700° fl + 0.710° fa R = 0.710° R - 0.700° H , 
0/ = 0.820° j + 0.580^, 0^ = -O.580? j7 + 0.820° 7 , 
= -0.650? + 0.760^, 0' = 0.760? + 0.650°,. (17) 

Thus, the lightest CP-even state 0# is the DM particle accounting for the cosmic ray anoma- 
lies. 0/, 0' and faj R are unstable but will contribute to the DM abundance in the freeze-in 
mechanism. The heavy metastable state is crucial to solve the small scale problem on 
power spectrum. Notice that we have arranged parameters to have < m-r, + m^ a , 
can not decay to 0' and <fin at two-body level. We emphasize that with suitable mass 0' 
might also constitute a component of DM today by forbidding its two body-decay to <pR 
and gravitino G. By the way, all the mixing factors are nearly democratic about 0.7, so for 
simplicity we may drop this factor in the following discussions, and the subscripts for 0j may 
be ignored since they will neither affect the discussions nor bring any misunderstanding. 

After the Z 3 symmetry breaking, the lightest Z 3 — odd state is unstable and decays to 
leptons through the heavy RHNs, which is described in Fig. [TJ The DM particles can decay 
via the operators Cijfafa^LiHy^^LjHy), and we are interested in the final sates containing \i 
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and r. At the leading order, such DM decays are described by the dimension-five operators 
C^(f)Vi{LjH u ) obtained from dimension- seven operators with one VEV for <pi and one VEV 
for H u . To show the close relation between DM decay and neutrino masses/TBM, we express 
the dimension-five operator coefficients into the light neutrino mass matrix elements. First, 
the Dirac neutrino mass matrix can be written as 

M LR = Y N (H° U ) = vsm(3 x (M,M,^ 3 ), (18) 

where Mi is the z-th column, and tan (3 = (H®) / (Hj) as in the MSSM, v = 
^({H°}) 2 + ((^)) 2 = 174 GeV. Using the seesaw formular M LL = M LR M^M^ R , we 
get the coefficients 



M x V M Nl J ijV smPM Nl 



Xx^ ( (M LL )ij\ ( 



wsin/3 V M x J \M J 



l N 2 



(19) 



where denotes v\ or V2, and the family universal couplings Xxi — 1 are assumed. This 
approximation is valid if M^ 2 dominates the seesaw contributions to Mn and M^ 2 ~ M^ 1 . 
Thus, the DM decays are closely related to the light neutrino mass matrix (elements), which 
will be studied in the next Section. We will show that the entries in the light neutrino mass 
matrix (Mll)^ are at the same order (about the heaviest neutrino mass) for i, j = 2, 3, 
while the other entries are much smaller (around the second heaviest neutrino mass or 
smaller) . 

The DM branch decay lifetime is 

r(<p R -> UiijHu) « 768tt 3 x 



(C|) 2 < 



3.6xlO-x(-^Yxf^^V M - 



10 15 GeVy V( M ^)u7 V^GeV 



5TeV\ 2 /2TeV 



x s, (20) 

H J V m 4>R J 

where we have taken tan /3 = 5 throughout this paper. The actual lifetime does not depend 
on it much since a larger tan/3 always gives sin /3 « 1. We keep \x<f> as an adjustable 
parameter to obtain the proper lifetime of (j) R , which will be chosen as 0.5 from then on. 
In order to generate the DM density via freeze- in mechanism, we choose M^jMx ~ 10 -11 . 
And then we explain the neutrino masses and mixings via the low-scale seesaw mechanism. 
Therefore, as pointed out in the Introduction, the crucial point to get such a long lifetime 
decaying DM is the combined factor Mn/Mx ~ 10~ 26 . 
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C. N2 Dominant Seesaw Mechanism and Cosmic Ray Anomalies 



Although our model can generate the suitable DM lifetime naturally, the dominant decays 
to the leptonic final states and the fittings of the PAMELA and Fermi-LAT data need further 
study. Especially, the decaying product should be dominated by the second and third families 
of charged leptons 2lJ. Note that the approximate \x — r symmetry is introduced to explain 
the light neutrino masses and mixings [22|, we suggest that the DM decay is related to the 
N 2 dominant seesaw mechanism which can explain neutrino TBM 12]. 

With approximate fi — r symmetry 
by four parameters 



22|, we obtain the general light Majorana mass matrix 



M 



LL 



m 



( X Y Y \ 

Y Z W 
\Y W Z ) 



(21) 



It predicts the maximal atmosphere mixing angle #23 = 7r/4 and #13 = 0, but leave the solar 
mixing angle #12 arbitrary. Taking sin 2 2# 12 = 8/9, the neutrino TBM is obtained 23J. The 
TBM Mil mass matrix only has three parameters since this fixed 612 is equivalent to a 
relation Z + W = X + Y . So we have 



fX Y 



M. 



LL 



m 



Y \ 



Y X+V Y-V 
\Y Y -V X + V ) 



(22) 



In the framework of seesaw mechanism with heavy RHN dominance, the crucial point of 
neutrino mixings is the specially aligned Dirac neutrino mass matrix (or say the Yukawa 
coupling matrix). Concretely speaking, the neutrino TBM can be understood by the aligned 
vacuum from an A4 discrete flavour symmetry breaking 24]. 

To explain why the DM decays dominant to muon and tau via neutrino physics, we 
modify the original Dirac Yukawa coupling ansatz used in Ref. [l^] as follows 



M 



LR 



(A 

A 

U 





-B 
B 



(23) 



To produce the realistic neutrino masses and mixings, we assume three RHNs with proper 
mass hierarchy < M^ 2 <C Mjv 3 so that the light neutrino mass matrix accommodates 
both the TBM and the \i + r dominated decay product of DM. Thus, the light neutrino 
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mass matrix is 



M; 



v 2 sin 2 /3 



MA/7 
Mjv, 



Mjv 2 



+ 



M^ 3 



m 



F F \ 
F X + F Y — V 
\Y Y -V X + V ) 



0{C 2 /M^ 



N 3/ 



(24) 



where X = Y = A 2 / (M^ttlq), and V = B 2 / (M/v 2 mo). The last term gives the subdominant 
contributions to Mll, but it is still important for the mass of the lightest neutrino. 

Now we show that the DM dominant decay channel to fi + r is a natural result for the N2 
dominant seesaw mechanism if the neutrino masses are normal hierarchy. Combining the 
DM decays with neutrino masses and TBM gives some contraints on the free parameters. 
First, it is obvious that Mll should be in the second RHN dominance. Next, with Eq. (EE 
we obtain three neutrino approximate masses 



2Vm 



2B 2 



M, 



N 2 



m U2 « 



3Xmo 



3A 2 



M, 



c 2 



(25) 



in a normal hierarchy form. Thus, the neutrino oscillation data Am 
2.40 x 10" 3 eV 2 suggest that 



and Am 2 } 



B 2 



A 2 



M N2 M 



5.4 : 1 



w 7.65 x 10" 5 eV 2 
(26) 



is valid when the N3 is sufficient heavy 25j. But from Eq. (fl9|) . the dimension- five operator 



coefficients are proportional to 1/Mjy . And then they disfavor large hierarchy Mjv 2 3> M^. 
So the hierarchy in Eq. ( 126]) is mainly due to the moderate relation A < B. Because B 2 /M^ 2 
will appear several times later, we fix it in the massless V\ limit (or say infinite m^ 3 limit) 



M, 



N 2 



Amfi/2 w 0.035 eV. 



(27) 



The DM 0^ decays to the SM fermions are the dominant primary source of comic ray 
since the lifetime of other states such as <fi is much short at the cosmic time scale. At tree 
level, the 0# three-body decay modes are 



j ■ 



iH u Vj , 



(28) 



and the corresponding lifetime estimation is given in Eq. (1201) . In Ref. 10], it has been 
explicitly simulated the electron spectra, and found that the spectra from direct hard leptons 
plus the soft contributions from cascade decays via the slpeptons and Higgs, Higgsinos are 
able to fit the PAMELA and Fermi-LAT experiments while not produce the anti-proton 
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excesses in the constrained MSSM. Moreover, for the two-body decay modes, <J)r decays to 
pure (two) neutrinos. The branch decay lifetime is approximately given by 



r((f) R -> UiVj 



57T X 



(Cf ? ) 2 m 3 



4>R 



m 4>R 
usin /? 



(29) 



which is about 20% of the one through three-body decays. Assuming a lifetime about 
5 x 10 26 s for three-body decay modes to explain the comic ray anomalies, we have r(0# — > 



10 s. The produced neutrino signals are potentially detectable with the upcoming 
IceCube neutrino observatory 26[, and the constraints on the DM models can be found in 
Ref. [22]. 



D. DM Density from Freeze-in Mechanism 

The decaying DM abundance can be produced naturally via the freeze-in mechanism (for 
a brief review, see Appendix |XJ) in our models. Let us explain the cosmological setup first 
since it is important to make freeze-in mechanism work. The initial relic density of <f>i should 
almost vanish, while (s)RHNs have the thermal density at least at T ~ M^. However, for 
the low-scale seesaw mechanism in the MSSM, the (s)RHNs only weakly interact with the 
plasma because the neutrino Dirac Yukawa couplings are small about 10~ 5 . Thus, after 
inflation the MSSM particles in the plasma alone cannot produce the thermal (s)RHNs. 
But this problem can be solved easily in the Next to the MSSM (NMSSM), where the 
superpotential term SNf can be introduced. Also, the suitable density for the (s)RHNs can 
be produced non-thermally by coupling them directly to the inflaton field. Thus, we assume 
the (s)RHNs in the plasma at the temperature T > M^.. But the initial densities for 0, are 
ignorable since they are SM singlets and only very weakly interacts with the (s)RHNs. 

However, during the decoupling of iVj, in the absence of inverse decay, the tiny branch 
decays or the scattering processes of Ni and Aj produces fa. To have the natural relic 
densities of fa via freeze-in mechanism, the typical couplings are required to be around 10~ 13 
for two-body decays and 10~ n for two to two scattering processes [ljj]. To be concrete, we 
give the relevant terms between (s)RHNs and (pi/ fa for freeze-in mechanism 

~ ~ ~ ~ ~ ~ 

-C DC^(fafa + fa<k)NiNi + —fafaN* + |i^| 2 + W^\ 2 + \F Nl \ 2 , 

-+CV?2 + v 2 ^ 1 )N l N l + C*M Ni (fxfaNiNt + /i.e.) + (...), (30) 

where we only consider the dominant terms, and dots denote many ignored terms, such as 
the supersymmetry breaking trilinear soft terms since the freeze-in amplitudes are controlled 
by M N . ^> Aij. Similarly, the scattering processes are also sub-dominated sinces they are 
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proportional to V{ which are several times smaller than Mn. in our model. In the precise 
calculations, one has to transform the interaction eigenstates to the mass eigenstates. For 
the DM state transformations, please see Eq. (fTTj) . In the following analysis we shall show 
that the DM relic density can indeed be obtained through the freeze-in mechanism, and the 
order-one mixing factor is not considered for simplicity. The mass eigenstates for jVj are 

iV i)± = -^(iV,±7V*), (31) 

where the squared mass eigenvalues are respectively given by M~ = M%. + m 2 ~ ± B^.. 
Here, m ~ are the soft mass square for the sRHNs. 

First, the freeze-in FIMPs from N it+ and iVj two-body decays N it+ — > and JVj — > Ni_<p 
as well as iVj. + — > Ni^tpm are in general kinetically forbidden. In fact, in the natural 
soft mass scale around (9(1 TeV), the mass splittings among N iy+ , „ and iVj are about 
M^usy/M^., which at most are tens of GeVs. Consequently there are no decay channels. 
In fact, it is required for proper DM relic density since the typical couplings given above 
are ~ M^./Mx 3> 10~ 13 . In short, these two-body decays must be forbidden (or at least 
suppressed sufficiently), otherwise, the freeze- in mechanism tends to over freeze DM(s) into 
the plasma. We have to point out that the above conclusion holds only when there are no 
mixings among the RHNs. If there exist the mixings in a complete model, we have to 
require that Max{ M N ., M N .}/M X ~ 1CT 13 since the mass splittings between the RHNs 
should be small enough to suppress the transition between JVj and Nj significantly. 

The scattering process JVjiVj — > (fiiR(fiiR from the second line of Eq. (I3"U|) can produce 
the phenomenologically important components ifim- The scattering process iVjiVj — > 4>4>ir 
can be studied similarly, so we will not present it in this paper. Numerically, the exact 
coincidence is a result of the dimension-five operators (pifaNiNj/Mx- Note that the RHNs 
have masses about 10 TeV, and Mx can be chosen a little bit smaller than Mqut, we obtain 
MfqjMx ~ 10~ n at the desired order. For the scattering processes, we consider the inter- 
action eigenstates as the mass eigenstates for simplicity since the mixings are democratic. 
The total cross sections are simply given by 

a(N l N l ^<P4 ) ) = ^\-V\l, x ~ ^~) . (32) 
oixs 

The function from phase space A(a, b, c) = (a — b — c) 2 — 46c (xi = mj/s) implies that only 
the lighter species (at least lighter than iVj) could be freezed into plasma with significant 
number densities. Numerically, the integral factor X[x, z] in Eq. (1A7j) is about 0.5, so the 
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relic densities of 4>m and are estimated to be at the right order 



6.0 x 

n ,h 2 Jl±[^±\\^ 



— &)(^)0(^) 2 . (33, 

where m$ denotes the mass of <f>m or 0. The relic density of gravitino G, which comes from 
the non-thermal production via the late decay, is about one order smaller than due 
to the mass ratio m^/mg ~ 10. We emphasize that this is not the final DM relic density, 
and the actual DM density is obtained by calculating all the processes with exact mixing 
factors. However, our results are enough to show that we can generate the correct DM relic 
density in our parameter space. 



III. PHENOMENOLOGICAL CONSEQUENCES 
A. Warmed G and Small Scale Problem 

In our model, both 0# and are generated with equal number densities via the freeze- 
in mechanism. Because is a relatively heavy metastable R p — odd state, it will decay 
and produce some relativistic particles. Thus, we can solve the small scale problem on 
power spectrum if the relativistic particle is the DM candidate like the LSP in the MSSM. 
If the comoving free-streaming scales of the relativistic particles, i.e., their motion in the 
comoving framework from their production time ij till to the matter and radiation equality 
era tEQ ~ 2.2 x 10 12 g can reach the small scale (9(0.1) Mpc, the power s pec trum on small 
scale can be reduced 28|. Such warm DM scenario was proposed in Ref. [l5[. 

To solve the small scale problem in our model, we require that have a proper mass 
(about 5 TeV in our example parameters) and a proper lifetime. Lifetime is fine in our 
model. Since is odd under Z 3 x R p , its leading decay mode is given by (as mentioned 
previously, in our interesting parameter space can not decay to 0' and 0r at two-body 
level) 

0->G + R ->... . (34) 

Because the decay rate is suppressed by 1/Mpj in gravity mediation, the lifetime can be 
sufficiently long about 10 — 1000 s, and can be even longer depending on mass splitting 
between the bosonic and fermionic states. However, if G is not the LSP, the above decay 
chain will produce a LSP such as neutralino. Thus, we have two viable solutions: (i) G 
itself is the LSP with mass about 0(100) GeV, then it is warmed enough to reduce power 
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spectrum on small scale; (B) G is not the LSP and has mass > (9(TeV). We require that 
G can produce the warm LSP via its decay while not forbid the two-body decay of 0. In 
short, these viable solutions do not conflicts with the parameter space in the MSSM. 

Because the late decay of G may spoil the successful predication of big bang primary 
nucleosynthesis (BBN), we consider G as the LSP for simplicity. In fact, this process can be 
regarded as a method of non-thermal production of G. The comoving free-streaming scale 
of a freely propagating particle can be calculated from the formular 15] 



R 



f 



ti 




a*Wi + W io« h/i + ,g (1 + ZeqY + igg + ZEQ) ) • < 35 > 

where and t^Q are the red shift and comic time at the matter-radiation equality era. 
Also, v is the current velocity of G 

Ti ttiq 

where To ~ 2.73 K, and Ej and Tj are respectively the energy and temperature when warm 
G is produced. According to Eq. ( )35|) . in order to explain the small-scale structure, v should 
take the value 1(T 8 - 10~ 7 |l5]. 

If G is light, vq is not dependent on its mass mg. Thus, with the proper mass for <fi in 
our example, we can indeed solve the small scale problem. Let us explain it in details. The 
two-body decay rate of <fi to its partner 0^ plus gravitino is calculated to be 29] 

1 m~ [ ( m \\ 
Tl = Wn^M\ l -\-^) j ' (3?) 

with the reduced Planck mass Mp = Mpi/i/87r ~ 2.4 x 10 18 GeV. Notice that this result is 
only valid when gravitino is much lighter than the mass splitting between particle and its 
superpartner, this is the exact situation needed in our model: G has a large velocity (warm 
enough) when it was produced. According to Eq. ( 1371) . the cosmological temperature (in the 

( -1/2 \ 1 / 2 

radiative dominant era) is given by Tj = 0.301#* ' Afpi/t/ , where g* is the effective 
relativistic degree of freedom in the plasma, and tj = 1/Tj ~ O(10)s with ~ 100 GeV. 
Then we can parametrize vq as follows 

/ 9* \W( ^ \ 3/2 f Am; \ 
V3.36/ \5TeVj \20TeV 2 J ' { } 

where Am| = m~ — m^ R , and we have used Ej = ^1 — m| fl /m~j ps As pointed 

out at the beginning, vq does not depend on mg explicitly, and the solution to the small-scale 
problem only depends on the mass of </>. 



v = 1.8 x 10~ 8 x 
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B. Baryon Asymmetry via Soft Leptogenesis 



Interestingly, we can also explain baryon asymmetry via the soft leptogensis, since the 
proper decay rate of 0r requires a low-scale seesaw mechanism with Mjq i ~ 10 4 GeV at 
least for i — 1, 2. Consequently, the new lepton number violation and CP violation in the 
supersymmetry breaking soft terms play major role in soft leptogensis. The soft leptogensis 
can work only well below < 10 £ 
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18j, which is the low bound on the right-handed 
neutrino mass in thermal leptogensis. And then reheating temperature is well below 10 9 
GeV as well. Thus, we can reduce the gravitino density produced by thermal scatterings 



in the thermal bath, and then the G late decay will not destroy BBN |30j. In short, the 
possible gravitino problem in the thermal leptogenesis can be solved. For a complete review, 
please see Ref. 3l|. In our model, the relevant Lagrangian for soft leptogenesis is 



1 



(m N2 + a* 2 )l\hI + ±Jf*_<y$y{M N% - a* 2 )l\hI 
i 



+ (Y£)*(N 2+ - N 2 _)L]Hl + -M 2 A 2 A 2 + h.c, 



(39) 



where Ai 2 = \Ai2\e l6A ^ ) and all the other soft terms have been taken real. Moreover, the 
SU(2) l gaugino mass M 2 is assumed to be real so that it will not induce large CP violation 
in the MSSM. 

In our model, the dominant contributions to the lepton number production come from 
the interference between the tree-level decays of iVj and the vertex corrections with gaugino 
running in the loop, which are given in Fig. |2j The original soft leptogensis relies on the 
self-energy corrections to JVj, and it is the small mass splitting (controlled by the bilinear 
soft terms B^.M^-Nf) between the two real degrees of freedom of A 7 , denoted with Ni± 
that resonantly enhances their CP-violation decays, see Fig. |2j To make the resonant effect 
large enough, < 1CT 3 Msusy must be fine-tuned to be very small 1?|, isj. Later, it was 
found that the vertex corrections to Ni decays with gaugino running in the loop contribute 
to the lepton number asymmetry in a very different way, and then the normal value of B^. 
is allowed [32j. This is important for our model to have successful soft leptogensis because 



its UV completion does not suppress B Ni . 

In the previous Section, we have considered the iV 2 dominant seesaw mechanism to pro- 
duce Mil. Corresponding to it, this dominance again dominantly generate the lepton asym- 
metry. This can be seen clearly from the explicit calculations of the lepton asymmetry 
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produced by a single Ni decays to lepton flavor a, using the procedure provided in Ref. 31] 

, _ 7 (U)+7(W-7(M-7(M 



A i2 \ 



3a 2 


m 1 


2 M 2 




\V N 


?M Ni 



log 



M| 



M| + Mjy. 



M 



A', 



sin ^ A BF (T), 



(40) 



where a 2 = g 2 /^ and g 2 is SU(2) L gauge coupling. This shows the N 2 dominant contribu- 
tions from Eqs. (I18|) and fl26|) . while the others are sub-dominant. Note that Y Nl ~ 10 -5 in 
our model, the contributions from self-energy corrections suppressed by \Y Ni \ 2 /a 2 are com- 
pletely ignorable. Abf(T), whose expression is given in Ref. 18| . denotes for the thermal 
effect in the thermal average of decay rates 7. Without it the above asymmetry vanishes 
at zero temperature field theory due to the exact cancellations between the fermionic and 
bosonic decay channels of iVj. By the way, our results are consistent with a previous work 
in Ref. 33j which used a different calculation method. 




< 

X 




FIG. 2: Lepton number and CP-violation decays of iV, with gaugino running in the vertex correction 
loop. Self-energy contributions are ignored since they are suppressed by the extra Yukawa couplings 
Y N \ 

The evolution of Ni lepton number violation decay and the evolution of a (a = e, /i, r) 
lepton flavor number are described by Boltzmann equations (BEs). Since the pure MSSM 
interaction conserves the charge A Q = (Bf + B s )/3 — (Lf + L s ) a , where = £ a + e a and 

L s is the total lepton number in scalar leptons. Similar definition applies to the Bf s . It is 



17 



convenient to study the evolution of density A a , and the coupled BEs are 36] 

z >y a (z) 



a 

Yk a = - e a (z)S a (z) + W a (z) (A* 



Y e J sH N 



(41) 
(42) 



where the derivative is on z = M^/T, and A^ = — Y*~ . In this crude set of BEs, we 
consider the AL = 1 two to two scattering processes that provide the CP asymmetry source, 
as well as the wash-out from the top quark and gauge boson interactions. Also, the flavor 
effect is kept for the low-scale soft leptogensis via the A matrix which expresses Yp a as the 
linear combination of Y Aa 34j, and via the C matrix which relates Yjj n to Y Aa 35| 



.4 



1 



207 



x 



/ 



V 



-64 
5 
5 



5 

-64 
5 



5 
5 



\ 



-64 / 



1^ 2 




(43) 



Our matrix is different from that in Refs. [33|, |36j since in this model the soft leptogenesis 
proceeds during the era T ~ ~ 10 4 GeV, where all the Yukawa couplings and the CKM 
mixings are in the chemical equilibrium. Also, the C matrix is included here since its entries 
are larger than the mixing entries in A matrix. Moreover, in the simplified BEs, the source 
term and wash-out terms can be rewritten analytically as follows 

KAz 



SJz] 



K 2 (z 



■K, a , WJz) 



(44) 



where the object lC a describes the degree of washout for a single flavor a 



/C, 



Hn m *MSSM 

\Y£\ 2 v 2 sm 2 f3/M N2 , 



(45) 



where m a is equal to \B\ 2 /M^ 2 for a = 2,3 while vanishes for a = 1. In the MSSM using 
g* = 228.75 at T ~ M N , one obtains m* MSSM w sin 2 x 1.58 x 10" 3 eV. Thus, /C 2 , 3 ~ 20, 
and then the soft leptogensis is in the strong wash out region 31_|, |34j . Finally, the sphaleron 
processes transform the survival lepton asymmetry into baryon asymmetry, eventually gives 
the baryon number density 

31 ^ a 



V MSSM 
1 B 



n B - rig 



(46) 



With the initial density of iVj in thermal equilibrium required by a successful freeze-in 
mechanism, we present the numerical solutions to the baryon asymmetry evolution in Fig. |3j 
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The observed baryon asymmetry Yb = (8.75 ± 0.23) x 10~ n [37| is generated with the 
following parameters: M N2 = 10 4 GeV, M 2 = 250 GeV and | jv 22 | = 300 GeV with phase 
9a 22 = — 1/4. As the RHN mass decreases, for instance, M^ 2 = 4 x 10 3 GeV, the baryon 
asymmetry tends to be overproduced. The reason is that the lepton asymmetry given in 
Eq. (140]) is proportional to 1/Mjy but linear to A^ 2 and M 2 . Because we have assumed that 
the LSP G has mass about 200 GeV, M 2 can not be too small. Therefore, we can choose a 
smaller |Ajv 22 | = 100 GeV with phase 6a 22 = ~ 1/5, °r we can fine-tune the phase of A N22 . 
Anyway, the observed baryon asymmetry can be obtained in the general parameter space. 



8.x 10" 
6.x 10" 
4.X10" 
2.x 10" 





10 12 14 



8 10 12 14 



FIG. 3: Baryon asymmetry Y AB (z) versus z = M N2 /T. Left: M N2 = 10 4 GeV, M 2 = 250 GeV and 
\A N22 \ = 300 GeV with 0a 22 = -1/4. Right: M N2 = 4 x 10 3 GeV, M 2 = 250 GeV, \A N22 \ = 100 
GeV, 6a 22 = —1/5. Initial density of N 2 is taken as thermal density. 



IV. DISCUSSIONS AND CONCLUSION 



Cosmic ray anomalies from the Fermi-LAT and PAMELA experiments can be naturally 
explained by the TeV-scale decaying DM with a very long lifetime ~ 10 26 s which decays 
dominantly to the muon and tau leptons. Note that the neutrino TBM can be realized 
elegantly via the \i — r symmetry, we conjectured that the DM decay is related to the 
neutrino physics. We considered the supersymmetric Standard Model with three right- 
handed neutrinos. To realize the decaying DM, we introduced a Z 3 discrete symmetry and 
two DM particles <f>\ and <p 2 . Because 4>i<p 2 can couple to the right-handed neutrinos via the 
dimension-five operators suppressed by the GUT scale Mqjjt, DM particle has a natural 
lifetime around r ~ 10 26 s if the seesaw scale is about 10 4 GeV. In particular, the DM 
particle will decay dominantly to the \x and r final states due to the N 2 dominant seesaw 
mechanism. Moreover, the DM relic density, which usually is a problem in decaying DM 
models, can be achieved naturally through the freeze-in mechanism with couplings typically 
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about 0(TeV)/MGUT- Simultaneously the small scale problem on the power spectrum can 
be solved since the metastable particles in the DM sector, which are also freezed into the 
thermal bath, can decay to the relativistic LSP in the supersymmetric SMs. Furthermore, 
we showed that the baryon asymmetry can be generated via the soft leptogenesis in a large 
region of the parameter space for supersymmetry breaking soft mass terms. 
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Appendix A: A Brief Review of Freeze-In Mechanism 

In this Appendix we will give a brief review of the freeze-in mechanism, but we shall 
formulate it differently. The basic idea of freeze-in mechanism is that in the BE for FIMP 
X, there is no inverse decay or scattering process to the mother particle that produces FIMP, 
as a result even small interaction rate is also successful in generating significant relic density 
for FIMP. As a starting point, the simple BE for X is (We consider the scattering process 
as an example, and the similar expression holds for decay.) 



where the Hubble constant at T = M A is Hi(T)\ Ma = 1.66y^T 2 /M p i\m a - I* 1 this paper, 

we use gP' s to denote the effective numbers of degree of freedom in the thermal bath at the 
freeze-in temperature T ~ M A , respectively for the entropy density s and energy density p. 

If A and B are assumed to be in thermal equilibrium during freeze-in, the BE reduces to 
the situation discussed in Ref. {l4j . The FIMP is produced dominantly at the temperature 
around the mass of heavier bath particles, when the bath particles still track their equilibrium 
distribution closely [ljj]. So the approximations are valid. In our paper, from Eq. (|4"T|) one 
finds that N and iV deviate from their equilibrium typically at Tf ~ (due to strong 

washout, sRHNs departure from equilibrium rather late), so the equilibrium approximation 
is also employed here. 

In our model, FIMP is freezed-in both from N± decays and their scatterings. First let 
us discus the decay. The yield of X is produced simply by integrating the right-hand side 
of Eq. flAlj) over z from to oo. A good property of freeze-in mechanism is that Y(x) 



(TL). 




sz^A + B -> X + C) 



Y A (z)Y B (z), 



(Al) 
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is insensitive to the UV physics, which is obvious from the integrand. So one can safely 
ignore the time of lower bound which may sensitive to inflation or reheating at the UV. For 
two-body decay, the thermally averaged decay rate is easily obtained analytically 



7 (4 X + C) = ^—^-z 2 Ki(z)T(A -> X + C), (A2) 







with g A the internal degrees of freedom of 4. Furthermore, let us reasonably assume that 
the freeze-in process lasts from z pa till z > 0{1) when generally the weakly interacting 
particle A decouples from thermal bath. Then one has 

r°° z 

Y(z> 10) « / dz—i{A -> X + C) 

135 qa M pl T(A^X + C) ^ 



8tt 3 (1.66)^V^ M i 
Then the relic density of X is given by 

n x h 2 ^ 4.50 x 10 25 x \ 2 —i^— — 



TF~i (A4) 

where we have typically used T(A — > X + C) = X 2 Ma/Stt and dropped the phase space 
factor. If multi thermal particles 4^ contribute to freeze-in, i should be summed over. Next 
we study the freeze-in mechanism through scattering processes. The thermally averaged 
scattering rate is formally given by 

7 (4 + B -» X + C) = 9A ff T 4 6 [ dxx 4 K 1 {x)X{l,x A ,x B ) 

lb7r J(m A +m B ) 2 /s 

x a(A + B -> X + C), (A5) 

where x = y/s/T, and then the integrand depends on z through its s dependence. Similarly, 
the final yield of X is obtained by integrating over z from to some large value 



r°° z 

Y x (oo) ^ / dz— - T (4 + B -> X + C) 
Jo sHx 

A5g A 9B M p i 

2 x 1.66 x 256^^^ 

i>oo i>oo 

X[x,z]= \ dz dxx 2 K 1 (x)E(x,z) ~ 0(1), (A6) 

where E(x,z) = (167rs)A(l, x A , xb)cf{x, z), and a(x,z) is the scattering cross section. Even- 
tually, the relic density is 

Q x h 2 Pa6.0 x 10 22 x -^L x ( ¥±) x X[x, z\. (A7) 

&<M \ m a) 
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